function Paras = Parameters()

    r=0.05;
    mu0=0.01;
    mu1=0.03;
    sigma=0.2;
    alpha=0;  %%%0.5%%%
    taue=0;  %%%0.4%%
    taui=0; %%%%0.05%%%
    I=20;
    p=18;
    x=1;
    pi0=0.5;
    eta=1;
    

    Paras = [ r mu0 mu1 sigma alpha taue taui I p x pi0 eta];

    
end